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Abstract. We study the influence of Stark broadening and stratification effects on Si I lines in the rapidly 
oscillating (roAp) star 10 Aql, where the Si I 6142.48 A and 6155.13 A lines are asymmetrical and shifted. First 
we have calculated Stark broadening parameters using the semiclassical perturbation method for three Si I lines: 
5950.2 A, 6142.48 A and 6155.13 A. We revised the synthetic spectrum calculation code taking into account both 
Stark width and shift for these lines. From the comparison of our calculations with the observations we found 
that Stark broadening + the stratification effect can explain asymmetry of the Si I 6142.48 A and 6155.13 A lines 
in the atmospere of roAp star 10 Aql. 
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1. Introduction 

The Stark broadening is the most significant pressure 
broadening mechanism for A and B stars and one has 
to take into account this effect in investigation, analysis 
and modeling of their atmospheres. In one of our previous 
works (Popovic et al. 2001) we have shown that the Stark 
effect may change the equivalent width of spectral lines 
by 10-45%, hence neglecting this mechanism, we may in- 
troduce a significant errors in abundance determinations. 
On the other hand, high resolution spectra allow us to 
study different broadening effects using line profiles. In the 
course of the abundance analyses of pecuhar (Ap) stars 
we noticed that some of the Sil lines are shifted rela- 
tive to the laboratory wavelength. Moreover, few strong 
lines mainly from the multiplets 3p^ ^D'^ — 5f^D and 
^Z?" — 5f^G have asymmetrical line profiles, in par- 
ticular the Si I 6155.13 A line. We found that this line is 
slightly shifted and asymmetrical even in the solar spec- 
trum, while in the hotter, e.g. Ap stars, the shift and asym- 
metry are more pronounced. 

The aim of this paper is to explain the asymmetry of 
Si I lines within the framework of the Stark broadening ef- 
fect. First we calculated the Stark broadening parameters 
for these lines, after that we calculated synthetic spectra 
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and compared them with the observations, and finally, we 
discuss the effects causing shift and asymmetry in the line 
profile of the Si 6155.13 and 6142.48 A lines. 



2. Observations 

For our analysis we used observations of one normal star 
HD 32115, two Ap stars HD 122970 and 10 Aql and 
Solar Flux Atlas by Kurucz et al. (1984). High resolution 
CCD spectra of 10 Aql and HD 122970 are described by 
Ryabchikova et al. (2000). High resolution CCD spectra 
(R«45000) of HD 32115 in the wavelength region 4000- 
9500 A were obtained with the coude-echell spectrome- 
ter mounted on the 2m 'Zeiss' telescope at Peak Terskol 
Observatory, Russia (see Bikmaev et al. 2002 for more 
details) . 

More Ap stars show peculiar line profiles of Si i lines 
but most stars have rather strong magnetic fields which 
distort line profiles through Zeeman splitting. Rather 
weak magnetic fields in Ap stars HD 122970 and 10 Aql 
allow us to ignore magnetic effects on line shape. 



3. The Stark broadening parameter calculation 

Calculations have been performed within the semiclas- 
sical perturbation formalism, developed and discussed 
in detail in Sahal-Brechot (1969ab). This formalism, as 
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well as the corresponding computer code, have been 
optimized and updated several times (see e.g. Sahal- 
Brechot, 1974; Dimitrijevic and Sahal- Brechot, 1984a, 
Dimitrijevic, 1996). 

Within this formalism, the full width of a neutral emit- 
ter isolated spectral line broadened by electron impacts 
can be expressed in terms of cross sections for elastic and 
inelastic processes as 



W, 



if 



2\jf 

2lTC 



and the corresponding line shift as 



A2 



2TTpdp sin 2(j)p 



(2) 



Here, \if is the wavelength of the line originating from the 
transition with the initial atomic energy level i and the fi- 
nal level /, c is the velocity of light, is the electron den- 
sity, f{v) is the Maxwellian velocity distribution function 
for electrons, m is the electron mass, k is the Boltzman 
constant, T is the temperature, and p denotes the impact 
parameter of the incoming electron. The inelastic cross 
section Ojji{v) is determined according to Chapter 3 in 
Sahal-Brechot (1969b), and elastic cross section (Jei ac- 
cording to Sahal-Brechot (1969a). The cut-offs, included 
in order to maintain for the unitarity of the iS-matrix, 
are described in Section 1 of Chapter 3 in Sahal-Brechot 
(1969a). 

The formulae for the ion-impact broadening param- 
eters are analogous to the formulae for electron-impact 
broadening. We note that the fact that the colliding ions 
would impact in the far wings should be checked, even for 
stellar densities. 



4. Line profile calculations 

Model atmosphere calculations as well as calculations of 
the absorption coefficients were made with the local ther- 
modynamical equilibrium (LTE) approximation. Model 
calculations were performed with the ATLAS9 code writ- 
ten by R.L. Kurucz (1993). The next step is the calculation 
of the outward flux at corresponding wavelengths points 
using the given model. For this we used the STARSP pro- 
gram written by V.V. Tsymbal (1996). In its current state 
this code includes the possibility of calculating a synthetic 
spectrum for an atmosphere with a vertical stratification 
of chemical elements. 

The computational scheme is as follows. For each line 
we find the central opacity as 



Tie -^n _h. 

gfife ^-(1-e '=^), 

mc p 



(3) 



where ai, is the mass absorption coefficient at fre- 
quency v, e is the electron charge, g is the statistical 
weight, fif is the oscillator strength for a given transition, 



X is the excitation energy, n is the number density of a cor- 
responding element in a given ionization stage multiplied 
by partition function, p is the density and h is the Planck 
constant. The last factor describes stimulated emission. 
The Doppler width is 



(4) 



In this expression niA is the mass of the absorber and 
is the microturbulent velocity. 

Next, we compute the total damping parameter 



7 = Irad + 1 Stark + Ineutral ■ 



(5) 



Here 7rad, 1 stark and ■jneutrai are the radiative. Stark 
and damping parameters due to neutral atom collisions re- 
spectively. The values oi ^j-ad, Ineutrai, cxcitation energy 
X and oscillator strength gf were taken from the Vienna 
Atomic Line Database (VALD) (Kupka et al. 1999). In the 
case of neutral atom broadening we assumed that perturb- 
ing particles are atoms of neutral hydrogen and helium 
only. This assumption is applicable to almost all types of 
stars due to high hydrogen and helium cosmic abundances. 
Usually this damping process is called Van der Waals 
broadening. The best theory for atomic hydrogen colli- 
sions which includes not only the Van der Waals potential 
is given in the papers by Anstee & O'Mara (1991) and 
Barklem & O'Mara (1998). Barklem, Piskunov & O'Mara 
(2000) provided \og^rieutrai calculations for about 5000 
lines of the neutral atoms and the first ions of many chem- 
ical elements. These damping parameters are included in 
VALD database per one perturbing particle and for a tem- 
perature of 10000 K. Unfortunately, no calculations exist 
for Si I lines. We shall discuss the competition between the 
broadenings caused by the Stark effect and neutral hydro- 
gen collisions in the atmospheres of our stars in Section 
5.2. 

In order to include Stark broadening effects we added 
the approximate formulas (see Eqs. (14) and (15) in Sec. 
5.1) in the code. The Stark width and shift are 



(e) , (p) , (Hell) 

1 Stark = Tstark^e + Tstark^P + Tstark ^Hell, 



(6) 



dstark = d^Stark^e + d^Stark^P + dftlrk ^Hell , (7) 

where rig. Hp and riHeii are the corresponding densities of 
electrons, protons and Hen ions respectively. The result- 
ing opacity profile is given by the Voigt function (Doppler 
-|- pressure broadening). 

Thus, at each point of a given spectral region (with res- 
olution 0.01 A for both lines) we computed line absorption 
coefficient as follows 



= a„V{u,a), 



(8) 



where is given by Eq. (4) and V{u, a) is a Voigt function 
with parameters 



V -1^0 + dstark 



(9) 



M. S. Dimitrijevic, T. Ryabchikova, L. C. Popovic, D. Shulyak, V. Tsymal: Stark broadening on Si I lines 



3 



7 



(10) 



The Stark shift dstark and the damping parameter 7 have 
been found from (2) and (11), respectively. 
The flux is given by the expression 



1 1 f°° 

H.ir.) = -2 S,E2{T,-t)dt+- J S,E2{t-T,)dt, (11 



where t„ is the optical depth, 5^ is the source function 
and En{x) is the exponential integral of the order n and 
argument x. The flux integral has been evaluated using 
matrix operators. 

5. Results 

5.1. Stark broadening data 

The atomic energy levels needed for Stark broadening 
calculations were taken from Martin and Zalubas (1983) 
and Moore (1971), but LS determination of hfD, hfG, 
Qs^P° and 7s^P° terms have been adopted according to 
Moore (1971) and therfore, in order to obtain a consis- 
tent set of data, energy levels from Moore (1971) have 
been used as the principal source. Oscillator strengths 
have been calculated by using the method of Bates 
and Damgaard (1949) and tables of Oertcl and Shomo 
(1968). For higher levels, the method described in van 
Regemorter, Binh Dy and Prud'homme (1979) has been 
apphed. 

The spectrum of neutral silicon is complex and not 
known well enough for a good calculation of the consid- 
ered lines. First of all upper and lower energy levels for 
6142.48 A and 6155.13 A lines are not known reliably. 
According to Striganov & Sventickij (1966) they belong 
to the 3p3 3£)0 „ 5 j3£, j^jjj 3jjo _ ^f3Q niultiplets re- 
spectively, while Moore ct al. (1966) stated that the lower 
level of the corresponding transitions is 3d^D°. It is also 
worth noting that in NIST (2002) the lower and upper 
levels of the corresponding transitions arc not specified. 
Consequently, we adopted the identification of Striganov 
& Sventitskij (1966) as the only one enabling the corre- 
sponding calculations. Moreover, the '3s3p^ "^D^ level is in 
fact a mixture of 39% of 3pnd^D° and 56% of 3s3p^ ^D^ 
states, and 5% is unknown. Also the 3pnd^D° state is a 
mixture of states with 3 < n < 12. In calculations, we as- 
sumed that we have 30% of 3p3d^D° state, and our checks 
show that the difference in final results is negligible if we 
assume that this state is involved up to 20% only. An addi- 
tional complication was that g levels, which might be very 
important for the perturbation of the considered 5/ levels, 
are unknown. In accordance with the decrease of distance 
between 5s, 5p, 5d and 5/ levels, we estimated that the 
distance from 5f^D or 5f^G term to the 5g levels should 
lie between 500 and 1500 cm~^. We checked results with- 
out 5g levels and with a Active 5g level at 500, 1000 and 
1500 cm~^ from the corresponding 5/ terms. In all cases, 
line widths differed by less than 1%, while the shift varies 
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Fig. 1. The analytic fit of He Il-impact broadening data 
for Si I 6142 A line: a) Stark width, and b) Stark shift. 
The shift and width are given for one perturber per cm^ 



within the limits of 3.5%. If the distance is 500 cm ^, the 
shift value e.g. at 10000 K for 6142.48 A line differs by 
about 20%. If we include a fictive 6g level 500 cm~^ dis- 
tant from 5g levels, the difference is negligible. Hence, in 
order to obtain the needed set of atomic energy levels, we 
adopted a fictive 5g level 1000 cm~^ distant from the 5/ 
levels and a fictive 6g level 500 cm~^ distant from the 5g 
level. Since the average estimated error of the semiclas- 
sical method is ±30%, due to additional approximations 
and uncertainties, we estimate the error bars of our results 
to be ±50%. 

Our results for electron-, proton-, and ionized helium- 
impact line widths and shifts for the three considered Si i 
spectral lines for a perturber density of 10^^ cm~'^ and 
temperatures T = 2,500 — 50,000 K, are shown in Table 
1. For perturber densities lower than those tabulated here, 
Stark broadening parameters vary linearly with perturber 
density. The nonlinear behaviour of Stark broadening pa- 
rameters at higher densities is the consequence of the in- 
fluence of Debye shielding and has been analyzed in detail 
in Dimitrijevic and Sahal-Brechot (1984b). 

In Table 1 the Stark widths and shifts for elect- 
ron-, proton- and Hell ion- impact broadening are pre- 
sented as a function of temperature for a density of 10^^ 



cm""*. After testing the density dependence of Stark pa- 
rameters, we have found that the width and shift are linear 
functions of density for perturber densities smaller than 
lO^^cm"^ and can be scaled by the simple formula: 



(W,dU = {W,d)o{^), 



(12) 



where {W,d)]\r are the width and shift at a perturber den- 
sity N (cm^^), and {W,d)o are width and shift given in 
Table 1, respectively. 
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Table 1. Stark broadening parameters for Sil spectral lines. This table shows electron-, proton-, and ionized helium- 
impact broadening parameters for Sil for a perturber density of 10^^ cm~^ and temperatures from 2,500 up to 50,000 
K. The quantity C (given in A cm~'^), when divided by the corresponding full width at half maximum, gives an 
estimate for the maximum perturber density for which tabulated data may be used. The asterisk identifies cases for 
which the collision volume multiplied by the perturber density (the condition for validity of the impact approximation) 
lies bectwen 0.1 and 0.5. For higher densities, the isolated line approximation used in the calculations breaks down. 
FWHM(A) denotes full lino width at half maximum m A, while SHIFT(A) denotes line shift in A. 



PERTURBERS ARE; 


ELECTRONS 


PROTONS 


HELIUM IONS 


TRANSITION 


T(K) 


FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


Sil 4s - 5p 


2500. 


0.509E-02 


0.333E-02 


0.160E-02 


0.701E-03 






ip-i D 


5000. 


0.581E-02 


0.409E-02 


0.172E-02 


0.898E-03 






5948.55 A 


10000. 


0.652E-02 


0.381E-02 


0.183E-02 


0.109E-02 


0.167E-02 


0.851E-03 


C= 0.36E+19 


20000. 


0.740E-02 


0.328E-02 


0.194E-02 


0.127E-02 


0.175E-02 


O.lOlE-02 




30000. 


0.811E-02 


0.271E-02 


0.202E-02 


0.139E-02 


0.181E-02 


O.lllE-02 




50000. 


0.903E-02 


0.212E-02 


0.213E-02 


0.153E-02 


0.188E-02 


0.123E-02 


Sil - 5/ 


2500. 


0.967E-01 


-0.576E-01 


0.175E-01 


-0.145E-01 


*0.145E-01 


-0.114E-01 




5000. 


0.108 


-0.624E-01 


0.196E-01 


-0.168E-01 


0.160E-01 


-0.134E-01 


6142.48 A 


10000. 


0.120 


-0.609E-01 


0.219E-01 


-0.192E-01 


0.178E-01 


-0.154E-01 


C= 0.52E-M6 


20000. 


0.132 


-0.472E-01 


0.247E-01 


-0.219E-01 


0.198E-01 


-0.176E-01 




30000. 


0.140 


-0.400E-01 


0.266E-01 


-0.235E-01 


0.211E-01 


-0.189E-01 




50000. 


0.147 


-0.326E-01 


0.294E-01 


-0.258E-01 


0.228E-01 


-0.208E-01 


Sil 3j3^ - 5/ 


2500. 


0.905E-01 


-0.622E-01 


0.184E-01 


-0.152E-01 


*0.151E-01 


-0.120E-01 




5000. 


0.102 


-0.707E-01 


0.204E-01 


-0.177E-01 


0.168E-01 


-0.141E-01 


6155.13 A 


10000. 


0.112 


-0.706E-01 


0.228E-01 


-0.203E-01 


0.187E-01 


-0.162E-01 


C= 0.20E+17 


20000. 


0.121 


-0.571E-01 


0.255E-01 


-0.231E-01 


0.208E-01 


-0.185E-01 




30000. 


0.129 


-0.482E-01 


0.273E-01 


-0.248E-01 


0.221E-01 


-0.200E-01 




50000. 


0.137 


-0.392E-01 


0.298E-01 


-0.272E-01 


0.240E-01 


-0.219E-01 




Fig. 2. A comparison between the observed Sil 6155 A 
line profile in the solar spectrum and synthetic spectra cal- 
culated with Stark widths and shifts from Table 1 (solid 
line) and with Stark widths calculated by the approxi- 
mate formulae and without Stark shift taken into account 
(dotted line). X- and Y-coordinates are wavelenghths and 
surface fluxes (normalized to unity). 

In order to simplify the use of Stark broadening data 
in the codes for stellar spectral synthesis, we have found 
an analytical expression for Stark widths and shifts 

^[k]=c,-{A + T% (13) 



Fig. 3. The same as in Fig. 2 for normal star HD 32115. 



— [A] =C2-(A + C-T^). (14) 



The constants Ci, C2, A, B and C are given in Table 2. We 
take T as T/{IQQQQK). The fits of our calculated data to 
the analytical functions given above are satisfactory (Fig. 
1) and the differences are smaller than the expected error 
of our calculations (ss ±50%) 
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Fig. 4. The same as in Fig. 2 for Ap star HD 122970. 
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0.70 LJ : : : : \ : : : : \ : : \ : : : : L 

6154.0 6154.5 6155.0 6155.5 6156.0 

Fig. 5. A comparison between the observed Sii 6155 A 
line profile in the spectrum of Ap star 10 Aql (thick line) 
and synthetic spectra calculated with Stark widths and 
shifts from Table 1 and Si abundance stratification (thin 
line), with the same Stark parameters but for homoge- 
neous Si distribution (dashed line), and with Stark width 
calculated by approximate formulae for the same strat- 
ification (dotted line). X- and Y-coordinates are wave- 
lenghths and surface fluxes (normalized to unity). 

5.2. The Stark broadening effect on the shape of Sii 
lines 

The asymmetry and shift of Sii 6142.48 A and (in par- 
ticular) 6155.13 A lines were observed, and it was obvious 
that the asymmetry is higher for the hotter stars. The 
most asymmetrical are A 6142.48 A and 6155.13 A, the 
strongest Unes of ^£,0 _ gjs^, ^^^^ 3^3 Zj^o _ ^fZQ 

multiplets, respectively. This line is slightly shifted and 
shows an asymmetry in the red wing even in a rather 
cool star like the Sun. It is not surprising if the Stark 
effect is responsible for the observed shifts and asymme- 
tries. From Table 1 one can see that neither widths nor 
shifts have a strong temperature dependence, but they de- 



1 .10 rr 




6141.0 6141.5 6142.0 6142.5 6143.0 



Fig. 6. The same as in Fig. 5 but for Sii 6142 A line. 



-3rT 
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10g(T501)o) 

Fig. 7. Si abundance distribution in the atmosphere of Ap 
star 10 Aql. 

pend linearly on the perturber density. If a line is strong 
enough so that line core and line wings are formed at dif- 
ferent atmospheric layers where electron density, differs 
significantly, then we should observe line asymmetry due 
to differnt shifts at different atmospheric layers even with 
a homogeneous vertical distribution of the absorbing ele- 
ment in stellar atmosphere (the case of the observed asym- 
metries in normal stars). Vertical abundance stratification 
makes this effect more pronounced, and it is observed in 
cool Ap stars, where an clement stratification takes place 
(Savanov et al. 2001, Ryabchikova et al. 2002, Wade et al. 
2001). 

Table 3 summarizes the model atmosphere parame- 
ters and rotational velocities for the stars of our sample. 
These data are taken from Ryabchikova et al. (2000) and 
Bikmaev et al. (2002). For the Sun corresponding param- 
eters were taken from Valenti & Piskunov (1996) except 
Teff. Si abundances used in the present calculations are 
given in the last column. 

In the Sun the most important broadening is due to 
collisions with neutral hydrogen and helium. In the hotter 



6 M. S. Dimitrijevic, T. Ryabchikova, L. C. Popovic, D. Shulyak, V. Tsymal: Stark broadening on Si I lines 

Table 3. The atmospheric parameters and rotational velocities of the investigated stars. 



Star name 


Tefl 


log g 




vsini 


log(Si/iVtot) 


or HD 






kms ^ 


kms~^ 




Sun 


5777 


4.44 


0.75 


1.6 


-4.49 


HD 122970 


6930 


4.11 


0.85 


5.0 


-4.45 


HD 32115 


7250 


4.20 


2.30 


9.0 


-4.65 


10 Aql 


7550 


4.00 


0.00 


5.0 


strat 




Table 2. The parameters A, B and C of the approximate 
formulae for Stark widths and shifts. 



6370.8 6371.0 6371.2 6371.4 6371.6 6371.8 6372.0 

Fig. 8. A comparison between the observed Sin 6371 A 
line profile in the spectrum of Ap star 10 Aql (double 
line) and the synthetic spectrum calculated with Si abun- 
dance distrubution shown in Fig. 7 (solid line) and Si ho- 
mogeneous abundance log(Si/A^tot) = -4.19 (dotted hue). 
X- and Y-coordinates are wavelcnghths and surface fluxes 
(normalized to unity), respectivelly 



stars the importance of the Stark broadening increases. P. 
Barklem (private communication) has provided us with 
the value of the broadening parameter due to collisions 
with atomic hydrogen for the Si i lines of our interest. This 
value per perturbing particle for a temperature 10000 K, 
logjneutrai =-6.63, is by 0.4 dcx higher than the corre- 
sponding value obtained from the Unsold (1955) approx- 
imation. According to P. Barklem his value may be over- 
estimated by 0.1 dcx which is proved by a comparison 
between the calculated and the observed blue wing of the 
solar Si I lines. The final value used in our calculations is 
log jneutrai =-6.75. As for the shift due to collisions with 
the atomic hydrogen, a corresponding theory is not well 
developped, and the value of the shift is rather uncertain 
(P. Barklem, private communication). The core and the 
wings of the strongest Si i A 6155 line are formed in the lay- 
ers with the temperatures 5000 and 5800 K respectively. In 
these layers the corresponding values for electron-, proton- 
and atomic hydrogen densities are 4.2 10^^ cm""^, 1.5 10^^ 
cm-3, 3.9 10^^ cm-3 - line core, and 2.5 10^^ cm-^, 9.3 
10^2 1.0 10^'^ cm-3 - line wings. The Stark width 

is 20 times smaller than ^ncuLrai at the depth of line core 
formation and 10 times smaller than ^neutral at the depth 
of line wing formation. In the solar atmosphere, the Stark 



Line 


5950.2 A 


6142.48 A 


6155.13 A 


WIDTH 


electrons 


Cl 


IE- 16 


IE- 14 


IE- 14 


A 


-0.33801648 


-0.87997007 


-0.88824105 


B 


0.12759724 


0.01704956 


0.01525209 


protons 


Cl 


lE-16 


lE-15 


lE-15 


A 


-0.81663954 


-0.77664655 


-0.76856840 


B 


0.01728187 


0.03913901 


0.03758103 


Hen 


Cl 


IE- 16 


IE- 15 


lE-15 


A 


-0.8:5337682 


-0.81987023 


-0.81114.334 


B 


0. 0129781;! 


0.02771059 


0.029 10282 


SHIFT 


electrons 


C2 


IE- 16 


lE-15 


lE-15 


A 


0.38028502 


-0.63997507 


-0.69208568 


B 


1.33496618 


0.88440353 


1.18033516 


C 


-0.02047206 


0.07950743 


0.04721782 


protons 


C2 


lE-16 


lE-15 


lE-15 


A 


-1.32564473 


1.12387002 


0.17032257 


B 


0.01917066 


0.02838789 


0.10512289 


C 


1.43394434 


-1.31786299 


-0.37318012 


Hen 


C2 


lE-16 


lE-15 


lE-15 


A 


-1.67370868 


0.14785217 


0.14840180 


B 


0.01490401 


0.10161093 


0.10439166 


c: 


1.7577;! 1«9 


-0.:!018:!:!:!9 


-0.:!1051296 



shift is 3 mA at the line core and about 20 mA for the 
line wings, which is not negligible. In the atmosphere of 
the hottest star of our sample, 10 Aql, the Stark width 
is of the order of ^neutral at the depth of line core for- 
mation and exceeds ^ineutrai by 5 times at the depth of 
line wing formation. Stark shifts are 10 and 100 mA, re- 
spectively. In 10 Aql the line core is formed at the layers 



with T=6500 K, iVc=1.4 10^^ cm'^, 7Vp=1.3 10^^ cm-^, 
while line wings are formed at the 
Np=1.2 10" 



7Vhi=5.9 lO^'"^ cm- 
layers with r=7500 K, N^=l.2 10^ 
cm-3, Nm=l.l 10^^ cm-^. 
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5.3. The Sun 

We first calculated Sil lines in the solar spectrum 
to check the Stark parameters. It was mentioned by 
Ryabchikova ct al. (2002) that available theoretical os- 
cillator strengths for the multiplets 3p^ — 5f^D and 
^D^ — 5f^G are overestimated by about 0.5 dex com- 
pared to the multiplet 4sip° - 5p^D (A 5948.55 A). The 
authors fitted the lines of - 5pD and 3p^ - 

5f^G multiplets using the solar spectrum. Corrected oscil- 
lator strengths - logg/=-1.42 (A 6142.48A) and log,g/=- 
0.77 (A 6155. IsA) were published in Bikmaev et al. (2002). 
The latter value was reduced to log(7/=-0.82 to get a bet- 
ter fit for the solar Sii 6155 line. We also reduced the 
oscillator strength (log5/=-2.25) of very weak Si line of 
the same multiplet, A 6155. 69A. We need to decrease both 
Stark widths and shifts for Sii 6142 and 6155 A lines by 
40% to get a reasonable fit of the synthetic line profiles to 
the solar ones. Due to relative weakness of the Sii 6142 
A line Stark effect is negligible but it is noticeable for 
the stronger line Sii 6155 A. Fig. 2 shows a comparison 
between the observed profile of the Sii 6155 A line in 
the solar spectrum and synthetic profiles calculated with 
(solid line ) and without (dashed line) taking into account 
the Stark width and in particular, shift. The strongest Sii 
5948.55 A line in our study is fully symmetric as expected 
taking into account an order of magnitude smaller values 
for both the Stark width and shift (see Table 1). 

With corrected Stark parameters we calculated Si i line 
profiles in spectra of other stars. 

5.4. HD 32115 

Spectra of this normal late A-type star have the smallest 
spectral resolution of our sample. As it may be seen from 
Fig. 3 the fit of the synthetic line profile of A 6155.13 A line 
to the observed one is much better when it was calculated 
with both Stark width and shift. As in the solar case, the 
Stark effect is negligible for the weaker Si i 6142.48 A line. 
To get a better fit of the line profiles we decrease the Si 
abundance by 0.13 dex compared to the results published 
by Bikmaev et al. (2002) which were based mainly on the 
equivalent width measurements. 

5.5. HD 122970 

This is one of the coolest Ap stars and is expected to 
have a chemically stratified atmosphere. Results of the 
abundance determinations show that except for the rare- 
earth elements (REE) abundance stratification is marginal 
if it exists at all (Ryabchikova et al. 2000). Spectral syn- 
thesis of Sii lines seems to support this conclusion. We 
could reproduce the asymmetric and shifted line profile of 
Sii 6155.13 A reasonably well, using the uniform distri- 
bution of Si and Stark broadening parameters presented 
here. Fig. 4 shows a comparison between the observed and 
computed line profiles for two cases similar to Fig. 3. The 



best fit of Sii line profiles was obtained with the same Si 
abundance as was derived by Ryabchikova et al. (2000). 

5.6. lOAqI 

10 Aql= HD 176232 is the hottest star in our sample. 
It has the most asymmetrical Sii 6155.13 A line profile, 
which could not be reproduced by any combination of 
Stark parameters in a homogeneous atmosphere (Fig. 5). 
The even weaker line, Sii 6142.48 A shows a noticeable 
line shift (Fig. 6). Ryabchikova et al. (2000) mentioned 
a possibility of Fe and the REE stratification in 10 Aql. 
Therefore we tried to find a simple distribution of Si in the 
atmosphere of 10 Aql by trial-and-error which would fit 
both Sii 6142.48 A and 6155.13 A lines. The simplest dis- 
tribution shown in Fig. 7 gives a reasonable fit to the ob- 
served profiles of both Sii lines (Figs. 5 and 6). Moreover, 
the same Si distribution seems to fit much better the pro- 
files of the strong Sin 6347, 6371 A spectral lines com- 
pared to the calculations with the homogeneous Si abun- 
dance (-4.19) obtained by Ryabchikova et al. (2000). It is 
beyond the scope of the present paper to derive an exact Si 
abundance profile in the atmosphere of 10 Aql. We would 
like to stress here that with the present Stark parameters 
a sensitivity of 6155.13 A line asymmetry to Si abundance 
changes in the stellar atmosphere can be successfully used 
in empirical studies of abundance stratification in the at- 
mospheres of cool Ap stars. 

6. Conclusions 

In order to discuss the contribution of the Stark broaden- 
ing effect to the asymmetry and the shift of Sii 6142.48 A 
and 6155.13 A lines we have calculated Stark broadening 
parameters for these lines by using the semiclasical per- 
turbation method. The obtained results are the first cal- 
culated data for the considered lines. In order to include 
Stark broadening data we have changed the STARSP code 
and we have synthesized the considered lines. From our 
analyzis we can conclude that: 

1) The Stark broadening effect is very important for 
these two lines. The contribution of electron impact is 
dominant but, impacts with protons and Hell ions should 
be taken into accout as well. 

2) The asymmetry as well as the shift of the Si I 6142.48 
A and 6155.13 A lines in many stars including the Sun can 
be explained by the Stark broadening effect. 

3) In hotter Ap stars, besides the Stark broadening ef- 
fect, the stratification plays a very important role in pro- 
ducing line asymmetry. The sensitivity of the line asym- 
metry to changes in the number of Si atoms through the 
stellar atmosphere can be used in abundance distribution 
studies. 
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